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Abstract. We present new, simultaneous, multifrequency observations of 45 out of the 55 candidate High Frequency 
Peakers (HFP) selected by Dallacasa et al. (2000), carried out 3 to 4 years after a first set of observations. 
Our sub-sample consists of 10 galaxies, 28 stellar objects ("quasars") and 7 unidentified sources. Both sets of 
observations are sensitive enough to allow the detection of variability at the 10% level or lower. While galaxies do 
not show significant variability, most quasars do. Seven of them no longer show the convex spectrum which is the 
defining property of Gigahertz Peaked Spectrum (GPS)/HFP sources and are interpreted as blazars caught by 
Dallacasa et al. (2000) during a flare, when a highly self-absorbed component dominated the emission. In general, 
the variability properties (amplitude, timescales, correlation between peak luminosity and peak frequency of 
the flaring component) of the quasar sub-sample resemble those of blazars. We thus conclude that most HFP 
candidates identified with quasars may well be flaring blazars. 
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1. Introduction 

The Compact Steep Spectrum (CSS) and Gigahertz 
Peaked Spectrum sources (GPS) are two classes of intrin- 
sically compact objects (linear size < 10-1 kpc) defined on 
' the basis of their spectral properties: the overall shape is 
■ convex with turnover frequencies between 100 MHz (CSS) 
and a few GHz (GPS) and the spectral index at high fre- 
quencies is steep (see O'Dea 1998 for a complete review). 

The currently accepted model (the youth scenario) re- 
lates the small linear size of GPS/CSS sources to their 
age, implying that they are the progenitors of extended ra- 
dio sources. Both kinematic studies (Owsianik & Conway 
1998; Owsianik, Conway & Polatidis 1998) and spectral 
analysis (Murgia et al. 1999) indicate ages of 10 3 -10 5 yr. 
On the other hand, various observations across the electro- 
magnetic spectrum (e.g. Fanti et al. 2000) failed to reveal 
the presence of a particularly dense medium that could 
prevent the expansion of the radio lobes in the intergalac- 
tic medium, as proposed in the frustration model (van 
Breugel 1984; Baum et al. 1990). So far no evidence of 
differences in the properties of the host galaxies between 
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GPS/CSS and FRII radio sources has been reported (e.g. 
see Fanti et al. 2000). 

Observational studies of the population of GPS/CSS 
sources have led to the discovery of an anti-correlation 
between the radio turnover frequency, v pi and the pro- 
jected angular size, (O'Dea & Baum 1997; Fanti et 
al. 2002). This relationship is expected if the peaked 
component is due to synchrotron self-absorption, where 
9 2 cx v p 5 ^ 2 , although free-free absorption could also play 
a role (Bicknell, Dopita & O'Dea 1997). In the youth sce- 
nario this correlation means that the youngest objects 
have the highest turnover frequencies and that the peak 
frequency is expected to move towards lower frequencies 
as the source expands and the energy density decreases. In 
principle, the highest the turnover frequency, the youngest 
the radio source is. 

The shape of the radio spectrum and the position of 
the turnover have been used as selection tools for this 
class of objects. Dallacasa et al. (2000) selected a sample 
of "bright" radio sources with turnover frequencies above 
5 GHz and called them High Frequency Peakers (HFPs). 
An early summary of their properties can be found in 
Dallacasa (2003). Simultaneous multifrequency VLA ob- 
servations were made to define the spectral shape of all 
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Fig. 1. Distribution of measured redshifts for the sample 
by Dallacasa et al. (2000). 

candidates. The comparison with literature data showed 
that flux density variability is not uncommon among the 
candidates and that also beamed radio sources like blazars 
can meet the selection criteria of the sample, e.g. when 
a flaring, strongly self-absorbed synchrotron component 
dominates the emission spectrum. A substantial contami- 
nation of the HFP sample by blazars may obviously lead 
astray analyses of HFP properties. 

In this paper we present the results of a second set 
of simultaneous multifrequency VLA observations of the 
sample of bright HFPs, aimed at discriminating between 
bona-fide HFPs and blazars, and at studying the spec- 
tral evolution of the sources. At 1.465 and 1.665 GHz, 
longer observations were carried out to search for extended 
emission around the point-like dominant component. It is 
known that a significant fraction (~ 10%) of GPS sources 
do have extended emission (e.g. Stanghellini 2003). Also 
polarization data have been taken and will be presented 
elsewhere. 

This paper is structured in the following way: Section 

2 describes the radio properties of the sample, the second- 
epoch VLA observations and the data reduction; Section 

3 presents the analysis of simultaneous radio spectra; 
Section 4 deals with the extended emission; in Section 
5 the variability properties are investigated; the main 
conclusions are summarized and discussed in Section 6. 
Throughout this paper, we adopt H = 50kms~ 1 Mpc~ 1 
and Q = 1 

2. The Bright HFP Sample 

2.1. Selection and radio properties of the sample 

The bright HFP sample of Dallacasa et al. (2000) was se- 
lected by cross-correlating the 87GB (Gregory et al. 1996) 



sources with S4.9GHZ > 300 mJy with the NVSS catalogue 
(Condon et al. 1998) at 1.4 GHz and picking out those 
with inverted spectra (a < —0.5, S oc v~ a ). Since the 
87GB and the NVSS surveys were carried out at differ- 
ent epochs, this sample is contaminated by flat-spectrum 
sources that happened to be flaring when they were ob- 
served at 4.9 GHz. The sample was then "cleaned" by 
means of simultaneous multifrequency VLA observations 
at 1.365, 1.665, 4.535, 4.985, 8.085, 8.485, 14.96 and 22.46 
GHz, leaving 55 sources whose single-epoch radio spec- 
trum peaks at frequencies ranging from a few GHz to 
about 22 GHz. These "first epoch" observations were car- 
ried out in 1998-1999. 

The final sample of HFP candidates comprises 11 
galaxies (including a type 1 Seyfert), 36 quasars, and 8 still 
unidentified sources (Dallacasa et al. 2002a). For compar- 
ison, galaxies and quasars are almost equally represented 
in GPS samples (Stanghellini et al. 1998) and most CSS 
sources are galaxies (Dallacasa 2002b, c). 36 objects out 
of 55 have known redshift and their distribution (Fig. ^) 
confirms the trend of quasars and galaxies to have barely 
overlapping redshift ranges. The redshift distribution of 
quasars shows a much slower decline with increasing red- 
shift at z > 2 than flat-spectrum quasars in the 2.7 GHz 
Parkes quarter Jansky sample (Jackson et al. 2002) with a 
similar flux-density limit (although at a lower frequency). 
This difference can be due, at least partly, to the more fa- 
vorable K-correction for high-z HFPs associated with their 
inverted spectra, and to the increasing variability am- 
plitude with increasing frequency (Impey & Neugebauer 
1988; Ciaramella et al. 2004), which also enhances the vis- 
ibility of distant quasars. 

2.2. Observations and data reduction 

In July 2002 we carried out multifrequency radio observa- 
tions of the 45 (out of 55) candidate HFPs that were visi- 
ble during the allocated observing time. We used the VLA 
in the B configuration, with a frequency coverage similar 
to that used for the sample definition. This sub-sample 
comprises 10 galaxies (including the type 1 Seyfert), 28 
stellar objects, henceforth referred to as "quasars" (in- 
cluding a BL Lac), and 7 unidentified sources. We used 
the intermediate frequencies (IFs) of 1.465 and 1.665 GHz 
in the L band, of 4.565 and 4.935 GHz in the C band, and 
of 8.085 and 8.465 GHz in the X band, while the standard 
VLA frequencies were used in the U band (14.96 GHz) and 
in the K band (22.46 GHz). The observing bandwidth was 
chosen to be 50 MHz per IF. Separate analysis for each IF 
in L, C and X bands was carried out in order to improve 
the spectral coverage of the data. 

Each source was observed typically for 1 min in C, 
X, U and K bands, and for 9 min in L band in a single 
snapshot, cycling through frequencies. The flux density 
measurements can thus be considered as simultaneous. 

Two scans were spent on the primary flux density 
calibrator 3C286, which was used also for the calibra- 
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Table 1. Integrated flux densities of candidates HFP sources at the 2002 epoch. 



J2000 


ID 


z 






Si. 4 


Si. 7 


S4.5 


5*4.9 


S 8 .i 


Ss.5 


5*15.0 


S22.5 


Name 






mjy 




mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


0003+2129 


G* 


0.4 


265 


c 


102 


120 


250 


253 


234 


227 


140 


86 


0005+0524 


Q 


1.887 


229 


c 


168 


186 


220 


213 


171 


166 


111 


82 


0037+0808 


G? 


>1.8 


292 


c 


101 


118 


281 


283 


267 


262 


190 


143 


0111+3906 


G 


0.668 


1324 


E 


476 


594 


1313 


1286 


978 


937 


507 


315 


0116+2422 






243 


c 


110 


124 


254 


258 


252 


248 




149 


0217+0144 


Q 


1.715 


1862 


E 


810 


802 


743 


745 


810 


810 


856 


838 


0329+3510 


Q* 




770 


E 


442 


462 


560 


555 


608 


613 


659 


702 


0357+2319 


Q* 




560 


c 


139 


133 


131 


130 


143 


144 


154 


168 


0428+3259 


G* 


0.3 


506 


E 


177 


202 


493 


506 


525 


514 


375 


263 


0519+0848 






278 


E 


279 


296 


446 


446 


435 


430 


420 


401 


0625+4440 


BL 




442 


c 


172 


184 


231 


231 


238 


238 


219 


210 


0638+5933 






591 


c 


277 


301 


590 


606 


668 


667 


620 


567 


0642+6758 


Q 


3.180 


474 


c 


245 


288 


429 


417 


331 


321 


203 


149 


0646+4451 


Q 


3.396 


1896 


c 


522 


612 


2580 


2770 


3709 


3757 


3691 


3318 


0650+6001 


Q 


0.455 


1236 


c 


517 


625 


1142 


1136 


989 


964 


671 


495 


1335+4542 


Q 


2.449 


735 


c 


294 


366 


797 


785 


613 


592 


359 


234 


1335+5844 






723 


c 


319 


396 


734 


725 


680 


671 


531 


253 


1407+2827 


G 


0.0769 


2362 


c 


910 


1100 


2469 


2463 


2114 


2050 


1139 


604 


1412+1334 






330 


c 


205 


236 


345 


337 


275 


267 


185 


130 


1424+2256 


Q 


3.626 


607 


c 


352 


414 


686 


669 


479 


460 


251 


145 


1430+1043 


Q 


1.710 


910 


c 


326 


402 


885 


882 


770 


752 


546 


385 


1457+0749 






241 


E 


360 


367 


262 


251 


194 


188 


141 


109 


1505+0326 


Q 


0.411 


929 


c 


413 


465 


740 


741 


721 


710 


665 


567 


1511+0518 


Syl 


0.084 


536 


c 


90 


115 


568 


607 


848 


861 


843 


617 


1526+6650 


Q 


3.02 


411 


c 


115 


145 


408 


410 


352 


341 


193 


107 


1603+1105 






270 


c? 


182 


188 


214 


214 


225 


225 


234 


217 


1616+0459 


Q 


3.197 


892 


c 


317 


382 


787 


771 


582 


559 


333 


212 


1623+6624 


G 


0.203 


298 


c 


159 


175 


295 


298 


287 


283 


224 


175 


1645+6330 


Q 


2.379 


513 


c 


303 


315 


510 


526 


615 


618 


596 


493 


1735+5049 


G? 




968 


c 


448 




925 


943 


934 


920 


740 


587 


1800+3848 


Q 


2.092 


791 


E 


313 


331 


702 


748 


1044 


1063 


1174 


1076 


1811+1704 


Q* 




691 


E 


545 


530 


494 


491 


508 


509 


499 


418 


1840+3900 


Q 


3.095 


203 


c 


125 


130 


169 


167 


161 


158 


134 


114 


1850+2825 


Q 


2.560 


1246 


c 


203 


223 


1096 


1185 


1540 


1541 


1318 


1045 


1855+3742 


G* 


0.5 


364 


c 


193 


191 


362 


345 


222 


212 


124 


91 


2021+0515 


Q* 




477 


c 


260 


302 


450 


441 


377 


368 


267 


191 


2024+1718 


s 


1.050? 


572 


E 


280 


304 


676 


704 


803 


800 


697 


569 


2101+0341 


Q 


1.013 


920 


E 


553 


610 


694 


678 


593 


583 


499 


478 


2123+0535 


Q 


1.878 


1896 


E 


1954 


2006 


2084 


2111 


2462 


2482 


2755 


2560 


2203+1007 


G* 


0.9 


319 


C 


124 


123 


312 


306 


243 


234 


129 


77 


2207+1652 


Q* 




551 


E 


223 




249 


247 


227 


223 


188 


163 


2212+2355 


S 




1182 


C 


601 




963 


979 


1033 


1028 


975 


915 


2257+0243 


Q 


2.081 


274 


C 


174 




274 


292 


439 


450 


558 


528 


2320+0513 


Q 


0.622 


1196 


E 


640 




657 


656 


720 


725 


806 


843 


2330+3348 


Q 


1.809 


558 


C 


370 


383 


407 


407 


457 


463 


532 


548 



* indicates the sources optically identified by Dallacasa et al. (2002a) . 



tion of the absolute orientation of the polarization vector. 
Secondary calibrators were observed for 1.5 min at each 
frequency about every 25 min; they were chosen aiming at 
minimizing the telescope slewing time. Accurate positions 
of the target sources were obtained from the JVAS cata- 
logue (Patnaik et al. 1992; Browne et al. 1998; Wilkinson 
et al. 1998). 

The data reduction followed the standard procedures 
for the VLA, implemented in the NRAO AIPS software. 



The L-band imaging has been quite complicated since a 
number of confusing sources fall within the primary beam, 
and an accurate flux density measurement could be ob- 
tained only once the confusing sources had been cleaned 
out. Generally, at least one iteration of phase-only self- 
calibration has been performed before the final imaging, 
although a few sources required several iterations. Some 
radio frequency interferences affected the 1.665 GHz data, 
so that a few sources could not be imaged and flux den- 
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sities could not be derived at this frequency. A Gaussian 
fit was performed on the final image by means of the task 
JMF1T, and the flux densities of the extended sources were 
determined with TVSTAT and 1MSTAT. 

The L band observations by Dallacasa et al. (2000) 
were not adequate to reveal possible extended emission 
given that were taken with the VLA in various configu- 
rations and were too short (1 min). The present deeper 
exposures detected extended emission for a considerable 
fraction of sources (see Sect. 4). 

The r.m.s. noise levels in the image plane are relevant 
only for flux densities of a few mJy and are generally con- 
sistent with the expected thermal noise. For our bright 
sources the main uncertainty comes from the amplitude 
calibration error, conservatively estimated to be (ler) 3% 
for the L, C and X bands, 5% for the U band, and 10% 
for the K band. 

The new multifrequency measurements are listed in 
columns 6-13 of Table □ together with the J2000 name 
(col. 1), the optical identification and the redshift taken 
from the NED database (columns 2 and 3, respectively), 
the 4.9 GHz flux density (col. 4) measured by Dallacasa et 
al. (2000); the labels C and E in col. 5 stand for compact 
or extended sources respectively, as seen in the L band 
images (see Sect. 4, Appendix). 

The time-lag between the the present data-set and the 
previous observations is of about 3 to 4 years. The years 
of observations of each source are specified in Table |21 

3. Spectral analysis 

In order to estimate the peak flux densities, S p , and fre- 
quencies, v p , of the sources, we have fitted the simultane- 
ous radio spectra at the two epochs using an hyperbolic 
function (Dallacasa et al. 2000) of the form: 

log(S) =a- [b 2 + {clogv - d) 2 ] 1 ' 2 (1) 

with the optically thin and thick spectral indices as 
asymptotes. The parameters of this function are related 
to S p and u p [log(f p ) = d/c, log(5 p ) = a — b], so that 
Eq. (JTJ can be rewritten as: 

log(S) = \og(S p ) + b-[b 2 + c 2 {\ogv - log vp) 2 } 1/2 . (2) 

The best fit values of v v and S p , obtained minimizing the 
chi-square function with the Minuit package (CERN li- 
braries), are reported, with their errors, in Table [3 The 
radio spectra of all sources are shown in Fig. [21 where the 
stars and the filled circles refer to the first and to the 
second epoch of simultaneous multifrequency VLA data, 
respectively, while the dashed and solid lines indicate the 
corresponding fits. 

A Kolmogorov-Smirnov (KS) test did not detect any 
significant difference among the distributions of the ob- 
served turnover frequencies (shown in Fig. |3J) and of the 
peak flux densities at the two epochs. This is not sur- 
prising since the time-lag is too short for an evolution of 
source spectra (see Sect. 5) to be detected. 



The observed peak frequency distributions of candi- 
date HFP galaxies and quasars are similar (left-hand panel 
of Fig. |J as a consequence of the selection criterion, but 
quasars are found up to much higher redshifts (Fig. HJ, 
implying much higher values of v p in the source frame 
(Fig. 01 right-hand panel). A similar behavior is found 
for GPS sources (Snellen 1997; Stanghellini et al. 1998), 
whose peak frequencies are, on average, about a factor of 
5 lower. 

4. Extended emission 

Within our sub-sample, 14 (31%) sources show some 
amount of extended emission (2 galaxies, 9 quasars and 
3 objects without optical identification) on scales ranging 
from 6 to 35 arcsec. In the Appendix we comment on such 
sources, based on our own observations and on published 
VLBA/VLA data. Our images are presented in Figs. Al- 
A10. The flux density of the extended emission has been 
measured as the difference between the integrated flux 
density over the whole source extension and the flux den- 
sity obtained using a point source model (see Table 0) . 
Weak emission below the surface brightness limit (~ 0.2 
mjy/beam) of our maps cannot be excluded. 

The total projected linear size of all these sources ex- 
ceeds the classical definition of a "compact" radio source, 
i.e. about 15-20 kpc (O'Dea 1998). The total flux density 
of the extended emission ranges between a few mJy and 
80 mJy, with 1.4 GHz luminosities ranging from ~ 10 24 
W/Hz (close to the upper end of the range for FR I 
sources) to ~ 10 27 W/Hz, well within the luminosity range 
of FR II sources. The unresolved to extended luminosity 
ratios are in the range 4 to 78. 

The fraction of sources with extended emission (~ 
30%) is larger than in GPS samples, where it is found 
to be ~ 10% (Stanghellini et al. 1990; Baum et al. 1990, 
Stanghellini et al., in preparation). However, no depen- 
dence of the fraction of sources with extended emission 
on peak frequency can be found within our sample: if we 
split it into two roughly equally populated sub-samples 
(with Up > 5.5 GHz and v p < 5.5, v p being the second 
epoch peak frequency in the observer's frame), we have 
approximately the same fraction of sources with extended 
emission in each. 

The presence of radio emission on scales of tens of 
kpc or even larger can indicate that we are really dealing 
with an evolved source but can also be reconciled with 
the youth scenario under the hypothesis of recurrent ac- 
tivity proposed by Baum et al. (1990), whereby a newly 
born source is propagating amidst the relic of previous 
large scale radio activity. At low frequencies (hundreds of 
MHz) the extended emission should be dominant, because 
of its steep spectrum. 

Seven out of the nine quasars with extended emis- 
sion have strongly variable spectra. Five of them did not 
show anymore a convex spectrum (see Sect. 5) when were 
re-observed, and therefore are not HFPs. On the other 
hand, the two galaxies with some amount of resolved emis- 
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Fig. 2. Radio spectra of sources [»!?„ (Jy) vs v (GHz)]. Stars and filled circles represent the first and the second epoch 
of simultaneous multifrequency VLA data respectively; the dashed and solid lines show the corresponding polynomial 
fits. 
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Fig. 2. (continued) 
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-. 1 — 1 1 1 1 1 1 1 — 1 1 1 1 1 1 1 — 1 1 1 1 II 1 — h= 


= 1 1 1 1 1 II 1 1 1 1 1 II 1 1 1 1 1 1 1 1 h 

I 2021 + 0515 Q* z=... 

r compact -_ 

1 — 1 1 1 1 II 1 — 1 1 1 1 II 1 — 1 1 1 1 1 1 1 — K 


: 2024+1718 S z=1.050 ! 
r extended ■= 

= — i — I 1 1 n 1 — I I 1 1 n 1 — I I 1 1 n 1 — h= 


2101 + 0341 Q z=1.013 
r extended - 

-. — i — I I 1 1 1 1 1 — I I 1 1 1 1 1 — I I 1 1 n 1 — h= 


: 2123 + 0535 Q z=l.B7B 

r extended - 

1 — 1 1 1 1 II 1 — 1 1 1 1 II 1 — 1 1 1 1 1 1 1 — K 


i 2203+1007 G* z = 0.9 

r compact -= 


2207+1652 Q* z = ... 
r extended - 


: 2212 + 2355 S z = ... 

r compact - 


; 1 1 1 1 II 1 1 1 1 1 II 1 1 1 1 1 II 1 H= 

: 2257 + 0243 Q z = 2.081 

r compact ■= 

1 — i — i — 1 — i — i — 1 — i — i — 1 — i — 


■ : 1 — MINI 1 — 1 1 1 III 1 — 1 1 1 1 II 1 — K 

2320 + 0513 Q z = 0.622 ! 
- extended - 

1 — i — i — 1 — i — i — 1 — i — i — 1 — i — 


= 1 — I I 1 1 II 1 — I I 1 1 II 1 — I I MM 1 — h 

: 2330 + 3348 Q z=1.809 '. 
r compact - 



1 10 100 1 10 100 1 10 100 



Fig. 



2. (continued) 
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Table 2. Peak frequencies and flux densities 



Name 






1998 


1999 


2002 


V 












Up 








Up 












TTlJv 


GHz 


TTlJv 


GHz 


mJv 


GHz 






0003+2129 


G* 


c 


272+5 


5.7+0.1 






257+5 


5.4+0.1 


3.02 


0.018 ± 0.009 


0005+0524 


Q 


c 






235+6 


4 13+0 09 


228+7 


3 40+0 09 


3.18 


0.170 ± 0.03 


0037+0808 


G? 


c 






287+5 


5.9+0.1 


287+6 


6.2+0.2 


0.46 


-0.05 ± 0.04 


0111+3906 


G 


E 






1333+28 


4.76+0.06 


1303+28 


4.68+0.07 


0.46 


0.01 ± 0.01 


0116+2422 




C 


245+4 


5.1 ±0.1 






266+9 


6.3+0.3 


5.59 




0217+0144 


Q 


E 






2557+105 


18+3 




flat 


238.27 




0329+3510 


Q* 


E 


768+13 


6.7+0.3 








flat 


21.60 




0357+2319 


Q* 


C 


637+15 


12+1 








flat 


420.12 




0428+3259 


G* 


E 


545+10 


7.3+0.2 






539+11 


6.8+0.2 


1.20 


0.02 ± 0.01 


0519+0848 




E 


>380 


>22 


>560 


>22 


448+7 


7.4+0.6 


62.80 




0625+4440 


BL 


C 


575+17 


13+2 






237+4 


7.4+1.0 


160.34 




0638+5933 




C 


700+17 


12+2 






675+13 


9.2+0.7 


0.83 




0642+6758 


Q 


c 


481+11 


4.5+0.1 






431+10 


4 08+0 08 


6.77 


0.10 ± 0.03 


0646+4451 


Q 


c 


3371+118 

*-J *J f A J L ±U 


15+2 






3982+92 


11.2+0.6 


36.91 


0.3 ± 0.1 


0650+6001 


Q 


c 


1990+19 


7.6+0.3 






1143+27 


5.2+0.2 


19.92 


0.11 ± 0.01 


1335+4542 


Q 


c 


737+18 


5.1+0.1 






793+17 


4.9+0.1 


1.23 


0.03 ± 0.02 


1335+5844 




c 


730+12 


6.0+0.2 






736+12 


5.5+0.1 


1.50 




1407+2827 


G 


c 


2320+36 


5 34+0 05 


2388+41 


5 20+0 09 


2484+49 


5 01+0 08 


1.77 


0.017+ 0.005 


1412+1334 




c 


337+7 


4.7+0.1 






345+8 


4.18+0.09 


1.03 




1424+2256 


Q 


c 


623+15 


4 13+0 07 






698+18 


3 94+0 06 


2.11 


0.05 ± 0.03 


1430+1043 


Q 


c 


905+14 


6.5+0.2 






887+39 


5.7+0.1 


4.19 


0.08 ± 0.02 


1457+0749 




E 






239+4 


5.4+0.3 


367+12 


1.7+0.3 


39.74 




1505+0326 


Q 


c 






937+14 


7.1+0.4 


744+11 


6.8+0.4 


19.65 


0.02 ± 0.04 


1511+0518 


Syl 


c 


778+18 


11.1+0.4 






903+20 

• 'y ''t - —V/ 


10.6+0.4 


8.28 


0.01 ± 0.01 


1526+6650 


Q 


c 


427+13 


5.7+0.1 






417+10 


5.5+0.1 


0.53 


0.04 ± 0.03 


1603+1105 




E 






277+4 


7.7+0.4 


229+9 


13+7 


16.65 




1616+0459 


Q 


c 






897+20 


4.7+0.1 


782+17 


4 63+0 09 


15.20 


0.02 ± 0.04 


1623+6624 


G 


c 


291+5 


6.0+0.2 






302+8 


6.0+0.2 


0.47 


0.00 ± 0.01 


1645+6330 


Q 


c 


628+18 


14+2 






629+17 


10.1+0.7 


2.74 


0.2 ± 0.1 


1735+5049 


G? 


c 


972±18 


6.4+0.2 






955+16 


6.3+0.3 


0.15 




1800+3848 


Q 


E 






1392+65 


17+3 


1226+49 


13+1 


1.91 


0.2 ± 0.2 


1811+1704 


Q* 


E 






809+19 


12+1 




flat 


58.74 




1840+3900 


Q 


c 






201+3 


5.7+0.5 


169+4 


5.2+0.4 


20.16 


0.1 ± 0.2 


1850+2825 


Q 


c 






1567+32 


9.1+0.3 


1591+32 

i • } ■ ' i — i — ' ) — 


9.5+0.3 


4.86 


-0.05+ 0.06 


1855+3742 


G* 


c 






426+42 


4.00+0.07 


423+310 


3.81+0.06 


1.32 


0.02 ± 0.01 


2021+0515 


Q* 


c 






497+13 


3 75+0 08 


446+10 


4.5+0.1 


5.92 




2024+1718 


S 


E 






845+32 


14+2 


807+16 


8.6+0.4 


9.58 


0.26 ± 0.06 


2101+0341 


Q 


E 






1635+90 


17+2 


737+212 


3.7+0.2 


119.87 


0.53 ± 0.02 


2123+0535 


Q 


E 






2392+103 


18+4 




flat 


27.49 




2203+1007 


G* 


C 






319+7 


4.86+0.07 


312+7 


5.0+0.1 


5.02 


-0.02+ 0.02 


2207+1652 


Q* 


E 






572+9 


7.4+0.3 


252+7 


3.5+0.3 


228.45 




2212+2355 


S 


C 






1379+13 


13+2 


1031+17 


9+1 


22.62 




2257+0243 


Q 


C 








>22 




>22 


1.85 




2320+0513 


Q 


E 






1176+20 


5.4+0.2 




flat 


86.24 




2330+3348 


Q 


C 






558+9 


5.6+0.3 




flat 


21.42 





sion on a kpc scale showed, at both epochs, the typical 
peaked spectrum. One is the well known J0111 + 3906 
(B0108 + 388) for which Baum et al. (1990) worked out 
the recurrent activity hypothesis. We are analyzing high 
resolution VLBA observations of J0428 + 3259 to study 
its pc scale morphology. 



5. Analysis of the sample 



The repeated high sensitivity simultaneous multifre- 
quency observations separated by 3 to 4 years allow us 
to investigate the variability properties of our candidate 
HFPs. This analysis may also help discriminating between 
truly young sources and flaring blazar components. 
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Name 


Sl.4 ( 


mJy) 


3l. T ( 


mJy) 


So band 


(mJy) 


<Sx band 


(mJy) 


LS (arcsec) 




peak 


ext 


peak 


ext 


peak 


ext 


peak 


ext 




J0111+3906 


477 


10 


594 


7 










17 


J0217+0144 


750 


75 


766 


59 


736 


20 


808 


8 


>6 (4.5 GHz) 


J0329+3510 


423 


80 


455 


65 


553 


5 






35 


J0428+3259 


167 


12 


198 


8 


492 


3 






>8 (4.5 GHz) 


J0519+0848 


268 


17 


287 


15 


444 


4 






10 


J1457+0749 


354 


13 


361 


12 










12 


J1603+1105 


170 


27 


175 


26 


213 


2 






94 


J1800+3848 


311 


4 


329 


5 










> 6 


J1811+1704 


537 


15 


523 


18 


488 


4 


494 


6 


22 


J2024+1718 


279 


12 


304 


10 










11 


J2101+0341 


550 


15 


609 


14 










22(EW)xl4(NS) 


J2123+0535 


1949 


14 


1999 


34 










11 


J2207+1652 


215 


50 






248 


12 


225 


6 


11 


J2320+0513 


623 


62 






656 


10 


721 


6 


19 (4.5 GHz) 



20 



i 1 i 1 i 1 r 



15 



epoch I 
. epoch II 



•Z, 10 



J , L 



12 16 20 



v v (GHz) 



Fig. 3. Histogram of observed turnover frequencies for the HFP sample at the two epochs. 



The convex form of the spectrum is the first property 
that should be preserved in the evolution of young sources. 
This does not happen for 7 sources (16% of the sample), 
all identified with quasars (25% of such objects), whose 
second epoch spectrum turns out to be flat. Such sources 
are labelled as "flat" in column 9 of Table |3 and are then 
classified as blazars. 

According to the physical model for the evolution of 
young radio sources by Begelman (1996), a newly gener- 
ated powerful jet injects energy into the ambient medium 
inflating a cocoon consisting of shocked jet material and 
shocked ambient matter. If the density of the surround- 
ing medium scales with radius, r, as p oc r _ ™, the linear 
size, LS, of the cocoon increases with the source age, t, as 
LS oc t e , with e = 3/(5 — n). The best determined rela- 
tionship for GPS/CSS sources relates the radio turnover 



frequency to the projected linear size v p oc (LS)~ S (O'Dea 
& Baum 1997). It follows that v p decreases with time as 
v p oc with A = 5 ■ e. The expected decrease of v v in 
a time interval At in the source frame (corresponding to 
AT = At(l + z) in the observer frame) is then: 



The slope of the observed correlation for CSS/GPS sources 
is 5 ~ 0.65 but O'Dea & Baum (1997) show that the 
evolution of v p with LS may be steeper than the statistical 
relation. Assuming Begelman's model, they found 6 ~ 1.2 
for n=2. 

The typical uncertainties associated with our estimates 
of Av p are at the several percent level. For n ranging from 
(uniform density) to 2 (isothermal distribution) and 
AT = 3-4 yr, a clear measurement of the peak frequency 
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Fig. 4. Observed (left) and rest frame (right) turnover frequency distributions for galaxies and quasars. 



decrease due to the source expansion could be achieved 
for source ages well below 100 yr. Taking into account that 
our sources were selected from the GB87 catalogue, based 
on observations carried out in 1986-1987, and therefore 
must have been older than 15-16 yr when our second 
set of observations were done, we conclude that sources 
with Az/p/ [z/ p Af(yr)] > 0.05 (last column of Table EJ are 
unlikely to be HFPs (because too young ages would be 
implied). We propose to classify them as blazars, whose 
observed evolution timescale are decreased by the Doppler 
factor. Note that Az^ p /[z/pAi(yr)] can be computed (and 
is reported in Table EJ) only for sources with measured 
redshift and measured peak frequencies at both epochs. 

In principle, an estimate of source ages could be de- 
rived from Eq. @. In practice, however, given the large 
uncertainties on Av p and the additional uncertainties on 
A (due to the poor knowledge of the parameter n char- 
acterizing the density profile of the medium), only rather 
uninteresting lower limits of a few tens of years can be 
set. On the other hand, interesting constraints would be 
provided by more accurate multifrequency measurements 
with a 2 or 3 times longer time-lag. 

Although the Kolmogorov-Smirnov (KS) test does not 
detect any significant difference in the distribution of 
turnover frequencies at the two epochs, it is clear from the 
last column of Table 2 that most sources have a smaller v v 
at the second epoch, consistent with an evolution towards 
lower optical depths. 

Flux density variability also helps distinguishing be- 
tween expanding young objects and aged beamed ob- 
jects (blazars) whose radio emission is dominated by 
a single knot in the jet. While HFP/GPS sources are 
thought to expand at mildly relativistic velocities, the 
large Doppler factors characterizing blazar jets boost the 
variability amplitudes and decrease the corresponding ob- 



served timescales. We have analyzed the variability of our 
sources in terms of the quantity 

v = j_f : (s, m -s,m\ (4) 

m of 

i=l 1 

where Si(i) and Su(i) are the flux densities at the i-th. 
frequency, measured in the first and second epoch, respec- 
tively, a is the error on Si(i)—Su(i), and m is the number 
of the sampled frequencies. 

We have noted that the second epoch flux densities at 
1.4 GHz are systematically higher by about 10% than the 
first epoch ones, while the flux densities at 1.7 GHz at the 
two epochs appear to match. This may be related to an 
unidentified problem during the 1.4 GHz observation of 
the primary calibrator 3C286. We have therefore chosen 
to exclude the 1.4 GHz measurements, so that in the above 
equation m = 7. 

For HFP sources, V depends essentially on only two 
free parameters, S p and v p (since the spectral slopes are 
not expected to vary), and should therefore have a \ 2 
distribution with 2 degrees of freedom. FigureJSJshows that 
indeed sources with low values of V obey such distribution. 
Since the probability that a source with V > 9 is extracted 
from such distribution is < 0.01, we conclude that sources 
with V > 9 are very likely blazars. On the other hand, 
blazars may well have small values of V, so that there is 
no guarantee that sources with V < 9 are all truly young. 
We caution however that these conclusions rely to some 
extent on our estimates of errors, dominated by calibration 
uncertainties: any significant over- or undcr-estimate of 
errors would impair the match of the distribution of V 
with the x 2 distribution with 2 degrees of freedom. 

As shown in Table [21 all 10 galaxies, including the 
Seyfert 1, have V < 9, while most quasars (17 out of 
28, including the 7 with flat second-epoch spectra) have 
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V > 9; of the 7 unidentified sources, 4 have V < 9. Thus, 
if we confine ourselves to sources with V < 9, we have an 
almost equal number of galaxies and quasars, similarly to 
what is found for GPS samples. The quasars with V < 9 
have higher median redshift (z me dian — 2.4) than those 
with V > 9, whose median redshift (z m edian — 1-7) is 
closer to that of flat spectrum radio quasars in the Parkes 
quarter- Jy sample (z modian ~ 1.4). 

An alternative approach may be to consider as blazar 
candidates all sources with variability significant at more 
than 3a. In this case the boundary between candidate 
HFP and likely blazars has to be set at V = 3 (as ap- 
propriate for a x 2 distribution with 1 degree of freedom) 
rather than at V = 9. If so, 9 more sources (3 galaxies, 
including the type 1 Seyfert, 5 quasars and 1 unidentified 
source) add to the list of candidate blazars. 

Another interesting issue is the relationship between 
peak frequencies and peak luminosities, L p . Opposite 
trends are predicted by current models for GPS sources: 
while Begelman's (1996) model implies a decrease of L p 
with decreasing i/ p , Snellen et al. (2000) predict an in- 
crease during the HFP/GPS phase of the source evolu- 
tion. In both cases, however, the dependence of L p on 
u p is weak, and therefore easily swamped by variances in 
both quantities. In the case of relativistic beaming with 
Doppler factor 6, the observed peak luminosities scale as 
S 3 , while the observed peak frequencies scale as 6, so that 
sources with equal intrinsic L p and v p , but different <5's are 
observed to have L p oc u p ; this relation is much steeper 
than the L p —u p relation predicted by both Begelman's and 
Snellen's models. 

We have investigated this relationship for all objects 
(see Fig. |5J, using first epoch data, when all sources but 
two had a well defined spectral peak. Although a clear 
correlation seems to be present, we must beware of the 
effect of the redshift distribution: as an example, the solid 
line in Fig. shows the distribution in the L p -v p plane of 
sources with equal observed S p = 300 mJy and v p — 5 GHz 
at different redshifts. Clearly the shape of the correla- 
tion is not far from that induced by the effect of redshift. 
Nevertheless, Kendall's partial correlation coefficient in- 
dicates a statistically significant positive correlation be- 
tween L p and v v after the influence of redshift has been 
eliminated. The probability of no correlation is 0.44% if 
we consider only stellar objects with V > 9 (which are 
probably blazars; we have 12 such objects with redshift, 
but for one of them the peak frequency is undefined) . For 
stellar objects with V < 9 the probability of no correlation 
increases to 6% so that the correlation is only marginally 
significant. For galaxies there is no indication of a signifi- 
cant correlation (probability of no correlation 38%). 

Models for GPS/HFP sources also predict a correlation 
between peak luminosity and linear size, whose slope is re- 
lated to the slope n of the density profile of the surround- 
ing medium. Unfortunately, our VLA observations provide 
only upper limits to the source sizes, since sources are un- 
resolved. A search of databases available in the VLBA 
Web site has yielded maps at various frequencies for 15 



out of the 18 sources with V < 9 and known redshift. Of 
these, 8 were resolved and for them we estimated the an- 
gular size. For the remaining 10 we adopted upper limits 
corresponding to the beam size. One source has not been 
included in the analysis because it has not a defined spec- 
tral peak at the first epoch of VLA observations. These 
data show no indication of a correlation between the linear 
size and the peak luminosity (probability of no correlation 
~ 0.7) and do not confirm, for our sample, the correlation 
between peak flux density and angular size (probability of 
no correlation ~ 0.6) reported by Snellen et al. (2000) for 
their GPS sample. 

6. Discussion and conclusions 

The spectral selection adopted by Dallacasa et al. (2000) 
is independent of optical identification and redshift. As 
a result, both galaxies and star-like objects ("quasars") 
are present in the original sample of bright HFP candi- 
dates. Quasars are almost three times more numerous than 
galaxies, consistent with the previously noted decrease of 
the galaxy to quasar ratio in GPS samples selected at 
increasing turnover frequency (Fanti et al. 1990; O'Dea 
1998; Stanghellini et al. 1998, 2003). 

However, our new multifrcquency observations, with a 
time lag of 3 to 4 years, indicate that the sample of can- 
didate HFP quasars is likely to be strongly contaminated 
by beamed objects, with variability properties consistent 
with those of blazars. 

Seven out of the 28 star-like sources (25%) no longer 
have a peaked spectrum, and 5 of them have extended 
emission. Such sources are therefore classified as blazars, 
caught by Dallacasa et al. (2000) during a flaring phase of 
a strongly self-absorbed component. 

Seven additional star-like sources show a variation of 
the peak frequency, significant at a > 3cr level, larger 
than expected if HFP/GPS sources expand at mildly rela- 
tivistic velocities (Av p /[v p At(yr)] > 0.05), suggesting the 
presence of relativistic beaming effects. The only galaxy 
(1407+2827) with a peak frequency variation significant 
at > 3a level has Ai/ p /[i/ p A<(yr)] = 0.017. 

Another criterion to discriminate between truly young 
sources and blazar is the variability index V [Eq. @] . The 
distribution of V for the less variable sources is close to 
the x 2 distribution with two degrees of freedom, as may 
be expected since variations of spectra of young sources 
are controlled by two parameters, S p and v v . If so, sources 
with V > 9 have a probability < 0.01 of being extracted 
from such a distribution and are therefore likely blazars. 
We have 20 sources with V > 9, 17 of which are star-like 
(including the 7 with flat spectrum at the second epoch 
and 4 of the 7 with Ais p /[i/ p At(yr)] > 0.05) and 3 unidenti- 
fied sources. The median r.m.s. variations for these sources 
are ~ 30% (27% for quasars only), close to those found for 
blazars (~ 32%; Ciaramella et al. 2004). 

The variability timescales in the frequency range 5-15 
GHz of bright blazars, derived from their structure func- 
tion, are ~ 2yr in the source frame (Hughes et al. 1992; 
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Fig. 5. Differential (left) and cumulative (right) distribution of V, compared with the corresponding \ 2 distributions 
for 2 degrees of freedom (dashed lines). The bins of the differential distribution correspond to A% 2 = 1. 




10 100 
u p (1+z) (GHz) 



Fig. 6. Peak luminosity versus peak frequency for all objects with measured redshift. Circles: quasars, squares: galaxies. 
Open symbols refer to sources with V < 9, filled symbols to sources with V > 9. The solid line shows the distribution 
in the L p -v p plane of sources with equal observed S p = 300 mJy and v v = 5 GHz at different rcdshifts. 



Ciaramella et al. 2004). To convert the time-lag between 
the two sets of observations into the mean time interval 
in the source frame, we need to divide it by (1 + z moan ). 
Only 11 of the 17 starlike sources with V > 9 have mea- 
sured redshift, with a median value of 1.7. We may thus 
expect that ~ 50% of blazars flaring at the moment of the 
first set of observations have recovered their baseline flat 
spectrum by the epoch when they where re-observed. For 
comparison, 7 out of 17 (i.e. ~ 41%) were indeed found to 
have a flat spectrum in the second observing run. 



There is evidence of a positive correlation between V 
and v p (probability of no correlation ~ 0.01), indicat- 
ing that the fraction of blazars is higher among objects 
with higher peak frequency. This is not surprising because 
young objects become increasingly rare since their life- 
times decrease with increasing v v . 

We have found a statistically significant positive cor- 
relation between the peak luminosity and rest-frame peak 
frequency of quasars, after having removed the effect of the 
redshift distribution (Kendall's partial correlation test). A 
relationship among these two quantities is expected both 
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in the case of GPS/HFP sources and of blazars, but with 
a much steeper slope for the latter sources. While previous 
studies enlightened a S p ~i>p correlation for a few individ- 
ual blazars (e.g., Stevens et al. 1996) our analysis indicates 
that it holds for flaring blazars as a class. 

There are thus various pieces of evidence, albeit cir- 
cumstantial, that most quasar HFP candidates are ac- 
tually objects where the effects of beaming are relevant 
(blazar like) even if this does not necessarily require that 
we are seeing evolved/old objects. 

While sources with V > 9 are likely blazars there is no 
guarantee that those with V < 9 are truly young sources. 
Some of them may really be blazars in a relatively qui- 
escent phase. The 3 such objects with A^ p /[^ p Ai(yr)] > 
0.05 may belong to this category. 

As for galaxies, all 10 HFP candidates have V < 9. 
Their peak luminosities and rest-frame peak frequencies 
do not show evidence of a correlation, not surprisingly 
given the smaller range spanned by both quantities (com- 
pared to the case of quasars) and the weaker dependence of 
L p on v p expected for GPS/HFP sources, easily swamped 
by intrinsic dispersions of both quantities. 

Two galaxies show extended emission compatible with 
recurrent activity. According to Begelman's model the 
synchrotron emission decreases with source age t as L oc 
t~ n with r\ = (n + 4)/ [4(5 — n)]. In the case of recur- 
rent activity, the HFP to extended emission age ratio 
can then be estimated as iiiFp/iext = (icxt/^HFp) 1 ^ 1 ', 
where L oxt and Lhfp are the total emitted radio lumi- 
nosities, that can be taken as proportional to the emitted 
flux densities at 1.4 GHz. In the case of the HFP com- 
ponent, the latter can be estimated extrapolating to 1.4 
GHz the peak flux density. The age ratio derived from 
Begelman's model depends on the slope n of the den- 
sity profile of the surrounding medium. Adopting n — 2 
we get, for J0111 + 3906, i 0Xt /< H FP = 9x 10 4 , and for 
J0428+3259 i xt/i H FP ~ 2 x 10 4 , consistent with the inter- 
mittency timescales suggested by Reynolds & Begelman 
(1997). Although any conclusion at this stage is prema- 
ture, this illustrates the potential of deeper low frequency 
observations to shed light on the nature of radio activity. 

In conclusion, while the variability properties of can- 
didate HFP galaxies may be consistent with expectations 
for truly young, possibly recurrent, sources, those of most 
quasars are consistent with those of blazars. If we set the 
boundary between HFP candidates and blazars at V = 3 
(see Sect. EJ, 3 galaxies would host a blazar nucleus, and 22 
out of the 28 star like sources and 4 out of the 7 unidentified 
sources would be blazar candidates. The median variabil- 
ity amplitude of sources with V > 3 is ~ 22%, slightly 
lower than, but still not far from that found for blazar 
samples. Milli-arcsec morphology at different frequencies, 
polarization measurements and possibly the study of the 
spectral index distributions in the optically thin regime 
can help distinguishing between flaring blazars and truly 
young objects. 

High resolution VLBA observations at two frequencies 
in the optically thin part of the spectrum (between 8.4 and 



43.2 GHz) of 50 objects of the original sample have already 
been acquired. The analysis of the data is in progress. 
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Appendix A: notes on the individual sources 

J0111+3906: well known galaxy at z = 0.668. We have 
detected extended emission on kpc scale (see Fig. lA.l|l . al- 
ready discovered and imaged by Baum et al. (1990), Taylor 
et al. (1996), Stanghcllini (2003), at different resolutions. 
The extended emission is in contrast with a recent origin 
for the radio activity in this source but can be explained 
if the radio source is recurrent (Baum et al. 1990). 
J0217+0144: quasar at z = 1.715. The source is slightly 
resolved in the NW-SE direction in the L band. At higher 
resolution, in the C and X bands, two opposite regions of 
emission with respect to the central compact component 
are detected (see Fig. IA.2|) . 

J0329+3510: stellar object. We see structure on arcsec- 
ond scale (lobes-hotspots). The extended emission is still 
detected, even if very weak, in the C band (see Fig. IA.3|) ; 
a jet like feature is emerging from the core towards the 
brightest lobe. 

J0428+3259: galaxy with estimated redshift z = 0.3 
(Dallacasa et al. 2002), based on the Hubble diagram 
in Snellen et al. (1996). The source is slightly resolved 
along the NW-SE axis at both 1.465 GHz and 1.665 GHz 
(Fig. IA.4|) . but no resolved emission is found at 5 GHz, 
due to sensitivity limitations. 

J0519+0848: optically unidentified. It is slightly resolved 
with a hint of weak emission in the NW direction (see 

Fig, rot. 

J1457+0749: optically unidentified. It is slightly resolved 
in the NW direction (see Fig. IA.6|) . 

J1603+1105: optically unidentified. The NE component 
can be either associated with the central object or can be 
an unrelated source (see Fig. IA.7|I ; according to Snellen 
(1998) there is a 6% chance of finding an unrelated NVSS 
radio source with a flux density of > 5 mJy within a radius 
of 100 arcsec from the target source. 
J1800+3848: quasar at z = 2.092. At 1.465 and 1.665 
GHz the source appears slightly resolved but the esti- 
mated flux density of the extended emission is of a few 
mJy only (see Fig. IA.8jl . 

J1811+1704: stellar object. A very weak component at 
NW is clearly visible both at 1.465 GHz and 1.665 GHz 
(see Fig. IA.8ll . 

J2024+1718: stellar object at z = 1.050. The emission 
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is extending from the core to the North for 11 arcsec (see 
Fig. rOjl . 

J2101+0341: quasar at z = 1.013. This source has a very 
complex structure. Two components are aligned with the 
NS axis for a total angular size of 14 arcsec. Furthermore 
there is a diffuse low-brightness emission in the NE direc- 
tion with angular size of 22 arcsec (see Fig. IA.9|l . 
J2123+0535: quasar at z = 1.878. There is a very weak 
emission toward the South of the strong compact compo- 
nent (see Fig. lA.lOj) . 

J2207+1652: stellar object. On arcsec scales it shows two 
components, the southern one being weaker and smaller. 
This component is still clearly visible in the X band (see 

Fig. En) . 

J2320+0513: quasar at z = 0.622. Extended emission is 
revealed in both the L and the C band. The maximum an- 
gular size in the NW-SE direction is of 19 arcsec. In the C 
band the extended emission has two irregular components 
on both sides of the unresolved nucleus, the southern one 
being weaker than the northern one (see Fig. IA.12|) . 



Patnaik, A.R., Browne, I.W.A., Wilkinson, P.N., et al. 1992, 

MNRAS, 254, 655 
Reynolds, C. S., & Begelman, M. C. 1997, ApJ, 487, L135 
Snellen, I.A.G., Bremer, M.N., Schilizzi, R.T. 1996, MNRAS, 

279, 1294 

Snellen, I.A.G. 1997, PhD Thesis, University of Leiden 
Snellen, I.A.G. , Schilizzi, R.T., & de Bruyn, A.G. 1998, A&AS, 
131, 435 

Snellen, I.A.G., Schilizzi, R.T., Miley, G.K., et al. 2000, 

MNRAS, 319, 445 
Stanghellini, C, Baum, S.A., O'Dea, CP., et al. 1990, A&A, 

233, 379 

Stanghellini, C, O'Dea, CP., Dallacasa, D., et al. 1998, A&AS, 
131, 303 

Stanghellini, C. 2003, PASA, 20, 118 

Stevens, J. A., Litchfield, S.J., Robson, E.I., et al. 1996, ApJ, 
466, 158 

Taylor, G.B., Readhead, A.C.S., & Pearson, T.J. 1996, ApJ, 
463, 95 

van Breugel, W., Heckman, T., & Miley, G. 1984, ApJ, 276, 79 
Wilkinson, P.N., Browne, I.W.A., Patnaik, A.R., et al. 1998, 
MNRAS, 300, 790 



References 

Baum, S. A., O'Dea, C P., de Bruyn, A. C, et al. 1990, A&A, 
232, 19 

Begelman, M.C 1996, in Proc. "Cygnus A - Study of radio 
Galaxy", Carilli C.L. & Harris D.E. eds., CUP, Cambridge, 
p. 209 

Bicknell, C, Dopita, M. A., & O'Dea, C P. 1997, ApJ, 485, 
112 

Browne I.W.A., Wilkinson P.N., Patnaik A.R., et al. 1998, 

MNRAS, 293, 257 
Ciaramella, A., Bongardo, C, Aller, H. D., et al. 2004, A&A, 

419, 485 

Condon J.J., Cotton W.D., Greisen E.W., et al. 1998, A.J. 115, 
1693 

Dallacasa, D., 2003 PASA, 20, 79 

Dallacasa, D., Falomo, R., & Stanghellini, C. 2002a, A&A, 382, 
53 

Dallacasa, D., Fanti, C, Giacintucci, S., et al. 2002b, A&A, 
389, 126 

Dallacasa, D., Stanghellini, C, Centonza, M., et al. 2000, A&A, 
363, 887 

Dallacasa, D., Tinti, S., Fanti, C, et al. 2002c, A&A, 389, 115 
Fanti, C, & Fanti, R. 2002,in Issues in Unification of 

Active Galactic Nuclei, ASP Conference Series 258, ed R. 

Maiolino, A. Marconi & N. Nagar (San Francisco: ASP), 

261 

Fanti, C, Pozzi, F., Fanti, R., et al. 2000, A&A, 358, 499 
Fanti, R., Fanti, C, Schilizzi, R.T., et al. 1990, A&A, 231, 333 
Gregory, P.O., Scott, W.K., Douglas, K., et al. 1996, ApJS, 
103, 427 

Hughes, P.A., Aller, H.D., & Aller, M.F. 1992, ApJ, 396, 469 
Impey, C. D., & Neugebauer, G. 1988, AJ, 95, 307 
Jackson, C.A., Wall, J.V., Shaver, P.A., et al. 2002, A&A, 386, 
97 

Murgia, M., Fanti, C, Fanti, R., et al. 1999, A&A, 345, 769 
O'Dea, CP., & Baum, S. A. 1997, AJ, 113, 148 
O'Dea, CP. 1998, PASP, 110, 493 
Owsianik, I., & Conway, J.E. 1998, A&A, 337, 69 
Owsianik, I., Conway, J.E., & Polatidis, A.G. 1998, A&A, 336, 
L37 



S. Tinti et al.: High Frequency Peakers: young radio sources or flaring blazars? 



15 



J01 1 1 +3906 - 1 464.900 MHZ 



J01 11+3906 - 1664.900 MHZ 



- 1 I I I I I I I 
f.c.= 0.15 mJy/b&am; PGak flux= 476.50 mJy/beam 


I I ~~ 






I I I I I I I 
f.c.= 0.12 mJy beam; Peak flux- 593.91 mJy/beam 


I I I 


- 

O 


o 




60 




— 


° 

— V O 






40 




* 


6 


















20 










o 

LU 












t/i 








o °WP 




o 

DC 





WW 


o _ 






< 














-20 




o 














— O 






-40 












-60 









I I - 









I 



100 SO 60 40 



20 -20 
ARC SEC 



100 80 60 40 



20 -20 -40 

ARC SEC 



Fig. A.l. J0111+3906; the restoring beam is in the Bottom Left Corner of each figure, the conturn levels are ±1, 2, 
4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096 times the first conturn (f.c.) corresponding to 3 times the r.m.s. noise 
of the image. 



J0217+0144- 1464.900 MHZ 



J0217+01 44 - 4935.1 00 MHZ 



f.c.= 0.22 mJy/beam; Peak flux- 749.75 mJy/beam 



60 - 
40 - 
20 - 
- 
-20 - 
-40 r 



© 



-50 
ARC SEC 



Fig. A.2. J0217+0144 



I 

- f.c.= 


I 

0.20 mJy/beam; Peak flux- 


I 

736.81 mJy/beam 










O 


<j 




















o 
















— f> 






















'J 










o 


O J 












o 






o 




El 


I 


i 









-10 -20 

ARC SEC 



16 



S. Tinti et al.: High Frequency Peakers: young radio sources or flaring blazars? 



J0329+3510 - 1464.900 MHZ 

_ f.c.= 0.1/ mjy/beam; Peak flux= 423.30 mjy/beam 



J0329+3510 - 4564.900 MHZ 



40 - 
20 - 
- 
■20 - 




-50 
ARC SEC 
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Fig. A.9. J2024+1718 (left) at 1.465 GHz and J2101+0341 (right) at 1.465 GHz 
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Fig. A.10. J2123+0535 
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Fig. A.ll. J2207+1652 
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